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Abstract. 

In order to analyze the dynamics of two languages in competition, one approach 
is to fit historical data on their numbers of speakers with a mathematical model in 
which the parameters are interpreted as the similarity between those languages and 
their relative status. Within this approach, we show here, on the basis of a detailed 
analysis and extensive calculations, the outcomes that can emerge for given values of 
these parameters. Contrary to previous results, it is possible that in the long term 
both languages coexist and survive. This happens only when there is a stable bilingual 
group, and this is possible only if the competing languages are sufficiently similar, in 
which case its occurrence is favoured by both similarity and status symmetry. 
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An aspect of globalization that alarms many is the replacement of local tongues by 
more hegemonic languages [1], a trend that has been investigated from multiple points 
of view, including that of Physics [21 [3l lU El El [7] . Where two languages compete, it 
is reasonable to assume that outcome is influenced by their relative status [6] (i.e. the 
perception by the speaker of the social and/or economic advantages each language offers) 
and their similarity [7] (or the complementary concept, interlinguistic distance [8]), but 
consideration of these factors has been hindered by difficulties in their operationalization. 

One of the earliest and simplest mathematical models of language shift was that of 
Abrams and Strogatz [6], who considered a stable population in which two languages 
with different statuses compete for speakers. This model, which involves the analysis 
of the evolution of the number of speakers along time, predicted that one of the 
languages would inevitably die out and was successfully fitted to historical data on 
competition between Scottish Gaelic and English, Welsh and English, and Quechua and 
Spanish, among other language pairs [6]. However, it did not take into account the 
possibility of bilingual individuals, a possibility that is of course realized in numerous 
multilingual societies. In Spain, for example, where Castilian Spanish is the official 
language throughout the State but in certain regions is co-official with another language 
(mainly Galician, Basque, Catalan or Valencian), individual bilingualism is common in 
communities with more than one co-official language. 

We recently showed that the historical evolution of the use of Galician and Castilian 
in Galicia (N.W. Spain) can be explained by a modified Abrams-Strogatz model that 
allows for bilingual as well as monolingual speakers of the competing languages, and 
which includes a parameter that represents the ease of bilingualism [7] (figure [1]). 
We considered a population in which monolingual X-speakers make up a fraction x, 
monolingual Y-speakers a fraction y, and bilinguals (B) a fraction b (with x + y + b = 1), 
and the dynamics of language change are accordingly described by the system 

yPyx + hPsx - x{PxY + Pxb), (la) 
xPxY + bPBY - yiPyx + Pyb), {Ih) 

xPxb + vPyb-KPby + Pbx)] (Ic) 

where Pxy denotes the probability of a monolingual X-speaker being replaced in the 
population by a monolingual Y-speaker, with analogous notation for the other possible 
replacements. The probability of a monolingual being replaced by a mono- or bilingual 
speaker of the other language is assumed to be proportional both to the status of the 
second language, i.e. the social and/or economic advantages it offers, and to a power of 
the proportion of the population that speak it. Thus, denoting by s the relative status 
of language X and by 1 — s that of language Y, 

PxB = c-k{l- s){l-x)\ (2a) 



dx 
dt 
&y_ 
&t 
db 
dt 



Pyb = c- ks{l - yY, 



{2b) 
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Figure 1. Flux among X-monolingual, Y-monolingual and bilingual groups. 

According to the model set forth in equations ([Ta|) - ([Tc|) and ((2a|) - ([2rf| . flux is governed 
both by the relative statuses of the competing languages (the relative social and/or 
economic advantages they offer) and by interlinguistic similarity fc, the probability 
that the disappearance of a monolingual speaker of one language will be compensated 
by the appearance of a bilingual rather than by a monolingual speaker of the other 
language. 



where c is a normalization factor related to the time scale, a is the power parameter, and 
k is the probability that the disappearance of a monolingual speaker of X (respectively 
Y) will be compensated by the appearance of a bilingual rather than by a monolingual 
speaker of Y (respectively X). We identify interlinguistic similarity with this parameter 



Note that since people do not actually forget their native tongue(s), the above 
model is a model of population renewal, rather than of individuals switching from one 
linguistic practice to another, except that it also allows monolinguals to become "born- 
again" bilinguals. When = it reduces to the Abrams-Strogatz model [6], or decays 
towards this model if b is initially nonzero. Note also that population growth does 
not invalidate the model, so long as the various linguistic groups are all affected in 
proportion to their size. 

Although the model sketched above adequately accounts for the Galician data up to 
the present [7], the question arises whether this situation is stable. More generally, what 
are the possible long-term outcomes of competition between two languages, and under 
what conditions do they come about? Namely, might two languages coexist stably? 







k. 
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To investigate these issues we have carried out extensive calculations, systematically 
varying s and k, to determine population states (x, y) that act as point attractors for 
the coupled system defined by equations ( fTal) -( |Tcl) and (l2al)-fl2dl) with x + y + b = 1; the 
problem is well-posed because, as it is easily shown, the velocity field on the boundary 
of the set of possible states (defined hj x > 0, y > 0, x + y < 1) never takes the 
system outside this set, whatever the values of s and k. K Initially, our calculations 
were performed as follows. For all {s,k) of the form (O.OSng, O.OSnfc) (0 < Ug^Uk < 20), 
10 000 starting states (xo,?/o) were randomly selected in the set of possible states, and 
the discretized form of the above system (fTaj) - (fTcj) was solved numerically for each, with 
a = 1.31, using a time step At = 0.01 (the value a = 1.31 was chosen because Abrams 
and Strogatz [6] found that this parameter was surprisingly constant, 1.31 ±0.25, when 
their model was fitted to 42 data sets concerning dissimilar language pairs). The 10 000 
calculations proceeded concurrently, and were halted as soon as all the state points had 
converged to within a circle of radius 10~^. Figures [2](a)-(d) show selected stages of 
this process for the case (s = 0.75, k = 0.3). However, it was soon found that there 
were values of (s, k) for which the state points did not converge, and examination of 
their distribution indicated that this was due to the existence of more than one point 
attractor, the attractor tended to by any particular state point depending on its starting 
state. In these cases the calculations were allowed to proceed for times several orders of 
magnitude longer than the average single-attractor convergence time, until all the state 
points were within two or three circles of radius 10~^ (figure [3] shows the result for the 
case s = 0.67, k = 0.77). 

Depending on s and k, the following five situations were observed to emerge: 

• I) There is just a single stable state at a; = 1 or ?/ = 1, i.e. one of the languages 
becomes extinct even though it may initially have been dominant. This is the 
behaviour illustrated in figures [2](a)-(d) and in figure Hl^a). 

• II) There are stable states at both x = 1 and y = 1; one of the languages dies out, 
but which one depends on the initial distribution (a;o,l/o), as illustrated in figure 

iKb). 

• III) There are stable states at a; = 1 and y = 1, together with a third that lies 
below the line x + y = 1 and thus corresponds to the presence of a stable bilingual 
group (figure m^c)). 

• IV) There is a stable state lying below the line x + y = 1 (i.e. with a non-empty 
bilingual group) together with just one stable monolingual state at a; = 1 (figures 
[la)-(d)) or y = 1 (figure SKd)). 

• V) There is just one stable state, and it includes a stable bilingual group (figure 
SKe)). 

Note that in cases II, III and IV there are sharp boundaries between the initial state 
zones leading to different outcomes; at these boundaries, an exogenous injection of just 
a few speakers into one group or another can determine whether a language lives or dies, 
as is illustrated in figures [5t^a) and[5](b). 
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Figure 2. Evolution of language dominance when languages X and Y 
compete. Example with s = 0.75 and k ~ 0.3. The graph axes represent the 
proportions of monolingual speakers of X and Y in the population, a, Random 
distribution of initial states, b-d, State distributions attained after respectively 100, 
300 and 1 500 steps of the computation. 

Figure Hf shows which of the above five situations was led to by each (s, k) value. 
The salient aspects of this map are that the stable existence of a bilingual group requires 
that k exceed a minimum value of about 0.35; that if /c > 0.6, then the less symmetric 
languages X and Y are statuswise, the larger k must be for stable bilingualism; and that 
when the two languages are moderately symmetric statuswise, which of them disappears 
depends on the initial sizes of the linguistic groups if they are essentially dissimilar 
{k < 0.4), but not if they are more similar (0.4 < k < 0.6). 

Returning to the case of Galician and Castilian in Galicia, the inclusion of recently 
published data |9] in the analysis (figure [5]^c)) corroborates our previous estimates [7] 
of Scaiician (0.26) and k (0.80), in fairly good agreement with lexicostatistics-based 
estimates of similarity among closely related Romance languages [10] ■ These values 
place Galicia in zone I of figure |l](f), and accordingly predict the eventual extinction 
of Galician. However, the close proximity of zone IV suggests, and figure [5](c) shows, 
that extinction is not imminent: in fact, it is predicted that by the end of the century 
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Figure 3. A case with two different attractors. State distributions attained 
after a, 200; b, 700; c, 6 000; and d, 10 000 computation steps when s = 0.67 
and k — 0.77 (the small regions shown in panels c and d contain all the points 
of the distributions). Depending on the initial proportions of X-monolinguals, Y- 
monolinguals and bilinguals, language Y dies out completely or survives, in the 
latter case mostly among bilinguals. Note the opening of a gap in the points of 
panel d, because the states evolve towards different attractors {x — l,y = 0) and 
{x = 0.6,2/ = 0.04). 



the population will be roughly equally divided between Castilian monolinguals and 
bilinguals. 

The model used in this work has two evident limitations. Firstly, the model does 
not consider possible alterations of the relative proportions of the linguistic groups due 
to immigration, emigration, or differential birth and/or death rates. Secondly, partially 
related to the previous one, is that the relative status of the two languages may well 
vary in time. Also, of course, it remains to be seen whether the definition of linguistic 
similarity in terms of the dynamics of the population renewal process corresponds to 
any purely linguistic concept; the assumption that similarity is symmetric, i.e. that 
-Pxs/[(1 — x)""] = PyB/[s(l —yY], is clearly an idealization; and there is the issue 

of whether two highly similar languages really have separate identities, especially when 
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Figure 4. The five possible stable situations determined by s and k. See 

text for descriptions of types I-V. Random distributions of initial states for typical 
cases of a, type I (s = 0.80, k = 0.65); b, type II (s = 0.40, k = 0.20); c, type III 
(s = 0.50, k = 0.65); d, type IV (s = 0.35, k = 0.75); and e, type V (s = 0.65, 
k — 0.85). Attractors are represented as large black spots. Where multiple attractors 
exist, initial states evolving towards different attractors are shown in different colours, 
f. The areas of the s, fc-plane occupied by systems of types I-V. The point k, s of the 
system Galician-Castilian is marked with a white spot. 
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Figure 5. Time evolution of single cases. Time dependence of the proportions 
of X-monolinguals (red), Y-monolinguals (green) and bilinguals (blue). Panels a and 
b illustrate for a type II system (s = 0.53, k = 0.55) how the fate of the system can 
depend critically on initial conditions, language Y becoming extinct if initially x = 0.84 
and y = 0.15 (a) but not ii x — 0.83 and y — 0.15 (b). c. Results of fitting the model 
to historical data for Galician (X, red), Castihan (Y, green) and bilinguals (B, blue) 
in Galicia (N.W. Spain). Although the fitted values of s and k, 0.26 and 0.80, place 
Galicia in zone I of figure |4ljf), thus predicting the eventual extinction of Galician, 
Galician-Castilian bilinguals are not expected to disappear within this century or the 
next. 
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they have similar status and one or both only survive among bilingual speakers. The 
present results nevertheless suggest that the competition between two languages does 
not inevitably lead to the extinction of one of them. 
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